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Interpretation of precision measurements of the cosmic microwave background (CMB) will require 
a detailed understanding of the recombination era, which determines such quantities as the acoustic 
oscillation scale and the Silk damping scale. This paper is the second in a series devoted to the 
subject of helium recombination, with a focus on two-photon processes in He I. The standard 
treatment of these processes includes only the spontaneous two-photon decay from the 2 1 5' level. 
We extend this treatment by including five additional effects, some of which have been suggested in 
recent papers but whose impact on He I recombination has not been fully quantified. These are: (i) 
stimulated two-photon decays; (ii) two-photon absorption of redshifted He I line radiation; (iii) two- 
photon decays from highly excited levels in He I (n 1 S and n 1 !), with n > 3); (iv) Raman scattering; 
and (v) the finite width of the 2 P° resonance. We find that effect (iii) is highly suppressed when 
one takes into account destructive interference between different intermediate states contributing to 
the two-photon decay amplitude. Overall, these effects are found to be insignificant: they modify 
the recombination history at the level of several parts in 10 4 . 

PACS numbers: 98.70.Vc, 95.30.Jx 



I. INTRODUCTION 

The anisotropy of the cosmic microwave background 
(CMB) has proven to be one of the most versatile and 
robust cosmological probes. The Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite has recently mea- 
sured these anisotropies at the percent level on de- 
gree scales [3, 0], and several experiments are ongoing 
or planned to make precise measurements of the po- 
larization and the sub-degree temperature fluctuations 
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CMB data, whose precision and robustness are so far 
unmatched by low-redshift observations, have provided 
some of the strongest tests of the standard cosmologi- 
cal model, including the adiabaticity and Gaussianity of 
the primordial perturbations and the spatial flatness of 
the universe. Most recently, the CMB has provided in- 
triguing evidence for departure of the spectrum of the 
primordial perturbations from scale invariance (n s < 1), 
as predicted by many models of inflation [ijj • 

The robustness of the CMB stems from the fact that 
the primary anisotropy can be calculated from first prin- 
ciples with reasonable computing time to a numerical ac- 
curacy of ~ 0.1% (i.e. good enough that this is not a lim- 
iting factor) [l8j]. The major exception to this statement 
is recombination, which affects the CMB anisotropy be- 
cause it determines the Thomson opacity and the visibil- 
ity function. The subject of cosmological recombination 
has a long history, with the early simple approximations 
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[TsL [2fj| being replaced by more sophistocated radiative 
transfer and multi- level atom analyses [ill. I22I I23L I24L I25I] . 
Several papers have appeared recently suggesting that 
the treatment of recombination in the current gener- 
ation of CMB anisotropy codes [H, HH] is incomplete 
[IE HE HE HE EE an d that the remaining errors may 
be large enough to be relevant for next-generation exper- 
iments such as Planck [HI]. It is especially worrisome that 
some of the errors in the standard recombination history, 
in particular helium recombination, are partially degen- 
erate with the scalar spectral index n^. a key parameter 
for constraining models of inflation [32, It is there- 
fore necessary to take a fresh look at the recombination 
problem. 

This paper ("Paper II") is the second in a series devoted 
to cosmological helium recombination. The first of these 
is Switzer & Hirata astro-ph/0702143, hereafter "Paper 
I," which re-examined helium recombination, taking into 
account the effects of semiforbidden and forbidden tran- 
sitions, spectral distortion feedback, and H I bound-free 
continuum opacity. We believe these are the major ef- 
fects in helium recombination that are not included in 
the standard treatment. This paper considers several re- 
visions to the standard treatment of two-photon transi- 
tions; these revisions do not have a major influence on 
helium recombination, but need to be included in order 
to establish that they are not important. The emphasis 
is on helium although some of the discussion (particu- 
larly that in Sees. HT1 and IIIip also applies to hydrogen. 
The third paper of the series (Switzer & Hirata astro- 
ph/0702145, hereafter "Paper III") will consider the ef- 
fects of 3 He scattering, electron scattering, rare decays, 
collsions, and peculiar velocities and summarize the ma- 
jor results. 
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The standard treatment of two-photon transitions in 
helium includes only the spontaneous two-photon decays 
from He I 2 1 S to the ground level l 1 ^, and their in- 
verse process, two-photon absorption. The first correc- 
tion considered in this paper is stimulated two-photon 
decay of 2 1 S to l 1 .!?, which was first analyzed by Chluba 
& Sunyaev [26] in the context of H I recombination. We 
reanalyze the effect here and also include two-photon ab- 
sorption of the spectral distortion as suggested by Kholu- 
penko & Ivanchik [3(3] , which delays recombination by re- 
exciting atoms. The stimulated two-photon transitions 
and absorption of the spectral distortion are found to 
play no significant role in He I or He II recombination, 
producing corrections to x e of the order of a few times 
1CT 5 . 

The second correction considered is the two-photon de- 
cay from highly excited levels (n^L, where n > 3 and 
L = 0,2), which was responsible for the largest correc- 
tion to recombination in the recent paper by Dubrovich 
& Grachev hereafter DG05. The treatment of such 
decays is a subtle issue because the two-photon spectrum 
from (for example) 3 1 D — > l 1 ^ contains a resonance as- 
sociated with the allowed sequence of one-photon tran- 
sitions 'i l D 2 1 P° -► l 1 ^. Since these "1+1" decays 
are already included in the level code one must be care- 
ful to distinguish which parts of the two-photon spec- 
trum should be added into the level code and which 
parts should be left out to avoid double-counting the rate. 
DG05 circumvented this difficulty by excluding the inter- 
mediate states associated with energetically allowed 1+1 
decays. This clearly avoids the double-counting prob- 
lem, but of course it affects the accuracy of the computed 
two-photon spectrum: in Sec. lIIII we will see that, partic- 
ularly for large n, DG05 overestimated the two-photon 
rate because they neglect destructive interference from 
the various intermediate states. 

In order to correctly implement two-photon rates from 
n > 3 in a level code, we must recall why they could be 
important even though they are much slower than the 
1+1 decays. The physical reason is that in a 1+1 de- 
cay, the higher photon is emitted in a He i n ,l P°-l 1 S 
line, and is likely to immediately re-excite another atom. 
There is no net production of the ground state He(l 1 S l ) 
except in the unlikely circumstance that the photon red- 
shifts out of the line or is absorbed by H I before it excites 
a He I atom. In contrast, the nonresonant two-photon 
decays in which neither photon is emitted within a He I 
line will produce a net gain of one ground state helium 
atom (except for the subtlety that one of the photons 
could later redshift into a He I line). Therefore, for the 
purposes of the level code, the way to distinguish "res- 
onant" (1+1) from "nonresonant" decays is not to make 
the distinction based on which intermediate state appears 
in the decay amplitude, but rather to impose a cutoff in 
frequency space: decays in which one of the photons is 
within Av cut of a He I n' 1 P°-l 1 S line are treated as res- 
onant (1+1), and the rest are nonresonant. The choice 
of Ai/ cu t (described in See.HV)) is arbitrary, reflecting the 



fact that the 1+1 decay is not a distinct physical process 
from two-photon decay - rather, the damping tails of the 
He in' 1 P^S line merge smoothly with the two-photon 
continua from all initial states that can decay to n' P°. 

Our approach to considering two-photon decays in this 
paper is to first consider the nonresonant decays for our 
choice of A^ cu t, and set an upper bound on how much 
they can speed up He I recombination by neglecting re- 
absorption of the nonresonant photons. The resonant 
two-photon decays (and the related processes of resonant 
two-photon absorption and resonant Raman scattering) 
can be considered as an alteration to the line profile of 
He I n^P^S, which is no longer well-described by a 
Voigt profile if one goes far enough out into the damping 
wings. Further, the 2 1 P°-1 1 S line now has a significant 
linewidth: for our choice of A^ cut , it requires ~ 0.02 Hub- 
ble times for a photon to redshift through the line (i.e. 
to redshift from frequency vn nc + A^ cut to Vu ne — Az/ cut ). 
Because of this, one must be careful about assuming the 
radiation field within the line is in steady state. All of 
these issues will be considered in Sec. El 

The outline in this paper is as follows. In Sec. Ull we 
consider the effect of stimulated two-photon decays from 
the n = 2 level (2 1 S) in He I and a related process, two- 
photon absorption of the spectral distortion. In Sec. IIII1 
we discuss the two-photon decay rates from highly ex- 
cited levels in He I (n > 3) and show that they were 
significantly overestimated by DG05. In order to evalu- 
ate the importance of the two-photon rates, we separate 
the two-photon spectrum into "nonresonant" and "reso- 
nant" pieces. The nonresonant contribution is considered 
in Sec. IIV} and the resonant contribution in Sec. El We 
conclude in Sec. IVII 

The notation in this paper is consistent with that in 
Paper I, but there are several new additions. Here we 
will denote the Rydberg constant by 1Z. The reduced 
matrix element of a spin k tensor operator is 
defined in accordance with Ref. [3J|. (This differs by a 
factor of i k from Ref. [35[, but is more convenient for our 
purposes because it makes the matrix elements real.) We 
will also use the symbol L> = max(i,l), which makes 
many appearances in our matrix elements. Spontaneous 
two-photon decay rates will be denoted by A, while the 
finite-temperature rates will be denoted T2 7 . Differential 
rates as a function of photon frequency or energy will be 
written dT2 1 /dv or dT2-y/dE. 



II. TWO-PHOTON DECAYS FROM n = 2 

The two-photon transitions from the metastable H I 
2s and He I 2 1 S levels are an important contribution 
to the recombination rates. It is usually assumed that 
stimulated two-photon emission plays a negligible role 
in the decay of the n = 2 states in H I and He I [25| . 
Chluba & Sunyaev [26] found that stimulated emission 
in H I 2s — > Is modifies TT and TE anisotropics at the 
percent level on small scales. We use a similar method to 
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include He I and He II stimulated two photon emission in 
addition to H I. We also include the re-absorption of the 
spectral distortion via two-photon excitation [3(j]. We do 
not find any significant effects in He I or He II. 

In this section, we present a general treatment of the 
two-photon decays from n = 2 levels, including stimu- 
lated emission and the effects of absorption of the spec- 
tral distortion. For He I n = 2 states, we consider first 
the usual two-photon decay: 

Hc(2 1 5') ^Hc(l 1 5 , ) + 7 spon + 7spon (1) 

and the two-photon excitation 

He(l 1 5)+ 7 bb + 7bb^He(2 1 5), (2) 

where 7 sp0 n refers to a spontaneously emitted photon and 
7bb refers to a photon drawn from the blackbody radia- 
tion. These two equations are the only ones considered in 
standard recombination codes, and they are typically in- 
cluded with a rate coefficient of An e i = 51 s" 1 for Eq. (Q]) 
and the detailed balance rate A H cie~ A - E/fcBTr for Eq. ([2]), 
where AE = E{2 1 S) - E^S). As pointed out for H I 
by Chluba & Sunyaev (2(|, one should also consider the 
analogous stimulated decays in He I: 



He(2 1 S') -» He(l 1 S') + 

7spon ~1~ Tstim 



and 



Re{2 1 S) -► He(l 1 5') 



' 7sti: 



' 7sti] 



(3) 



(4) 



where 7 s tim means that the photon's emission is stimu- 
lated. [Chluba & Sunyaev [26| replaced Ahi in their level 
code with the sum of rates for Eqs. {J), |(3|), and Q.] 
Note that it is not self-consistent to leave out these re- 
actions, since the reverse reaction of Eq. (f2|) is not just 
Eq. (fl]) but rather the combination of Eqs. (H|), ©, and 

©• 

If there is a spectral distortion from redshifted He I line 
photons, one should also consider the possibility of two- 
photon absorption of a thermal photon and a distortion 
photon 



He(l 1 5) + 7 bb + 7dist - He(2 1 5), 



(5) 



where 7dist refers to a spectral distortion photon. In prin- 
ciple, there is an additional contribution where both ab- 
sorbed photons come from the spectral distortion. This is 
negligible since the blackbody spectrum dominates over 
the spectral distortion for photons with hu < AE(2 1 S)/2 
except at z < 1400 when He I recombination is finished 
(xHeii < 10 -14 ). All of these equations have analogues 
in H I and He II. 

The two-photon decay rate is 



%2 1 S^1 1 S 



= A 



Hcl 



I/ llS-2ls/ 2 



(C) 



where < 
to 



{v) is the two-photon emission profile, normalized 



'llS-2ls/ 2 



dv = 1, 



(7) 



and the frequency of the higher-frequency photon is v' = 
vi 1 s-2 1 s~v- Note that decay term [1 +hf(v)} [1 +N{v')\ 
can be expanded to give a spontaneous piece, a singly 
stimulated piece N{v) +N{v'), and a doubly stimulated 
piece J\f(v)J\f(v'). The phase space density for the higher- 
energy photon is much less than unity: at z = 2600 it is 
2 x 10 -4 at the midpoint of the H I spectrum, fi s _2 S /2, 
and it is even less above the midpoint, for He I, or for 
lower redshifts). Therefore we make the replacement in 
the downward rate 1 + N{v') — > 1, i.e. we neglect stim- 
ulated emission of the higher-energy photon. Similarly 
since the spectral distortion phase space density is <C 1, 
we may replace 1 +J\f(v) — > 1 + A/t>b(^)- This enables us 
to write Eq. ([6]) as 



. (thermal) 



• (nonthermal) 
C 2 1 S->l 1 S ' 



(8) 



where 



. (thermal) 
X 2 1 S^1 1 S 



A 



Hcl 



Jy llS-2ls/ 2 



4>{v) dv 



x a^S - X!i S e 



-hv 1 i s _ 2 i s /kT r 



(9) 



and 



. (nonthermal) 



-A 



chvlkT r _ I 



XA/clistlV^l 1 ,! 



(10) 



[Here A/dist(^') is the distortion phase space density de- 
fined by taking the actual phase space density and sub- 
tracting the blackbody contribution.] 

In the level code, the profiles for H I and He n are based 
on the fits by Nussbaumer and Schmutz [3(|, and for He I 
we use the fit to Drake described in Appendix A of 
Paper I. 

The results of including Eqs. Q and (fT0|) in the level 
code are shown in Fig. [TJ We can see that for He I and 
He II the effect is very small - only a few times 10~ 5 . A 
larger effect in Ax e during hydrogen recombination was 
found by [2£|, for several reasons. First, the absolute 
abundance of hydrogen is greater, so a similar fractional 
change in its recombination history leads to a larger ef- 
fect. Second, the 2s and 2p levels in H I are essentially 
degenerate, whereas in He I the 2 1 P° level lies 0.6 eV 
above the 2 1 S level. This changes the shape of the two- 
photon spectrum at low frequencies, where 4>{v) oc v 3 in 
He I (due purely to the available phase space for emit- 
ting a low energy photon) as opposed to <fi(v) oc v in 
H I (where there is a pole corresponding to the 2p inter- 
mediate state in the matrix element at zero frequency). 
Stimulated emission and re-absorption of the spectral dis- 
tortion will play a larger role in the case of H I where the 
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Ax e from modifications of two-photon transitions from n=2 

4x1 0" 5 i 1 1 1 1 1 1 




2000 3000 4000 5000 6000 



FIG. 1: A comparison of the effect of stimulated two-photon 
emission and nonthermal two-photon absorption relative to 
the reference model. The net effect is a delay in recombi- 
nation, Ai e > 0. The two peaks correspond to the effect 
on He ii recombination (z ~ 5500) and He I recombination 
(z ~ 2500). Note that in both cases the effect on the recom- 
bination history is small, i.e. a few parts in 10 5 . 



two-photon spectrum has more probability at the ends 
of the spectrum. A third reason is that, due largely 
to the lower abundance of helium versus hydrogen, the 
He I 2 1 P°-1 1 S optical depth during He I recombination 
is much less than that of Lya during H I recombination; 
therefore the importance of two-photon decays relative 
to resonance escape is less for He I than for H I. 



III. TWO-PHOTON RATES FROM n > 3 

In this section, we consider the effect of two-photon de- 
cays from the higher excited levels of He I (n > 3) . Such 
decays were discussed by DG05 as a potential means to 
dramatically speed up He I recombination. The inclusion 
of these decays introduces a new subtlety, however, which 
is not present for the n = 2 two-photon decay 2 1 5 — > IS. 
The 2 1 S level does not have any allowed decay routes, so 
it is correct to take a multilevel atom code and add in a 
new rate ^s-^^s for this decay. In contrast, the higher 
levels in He I (n 1 L, n > 3, L <G {0,2}) that can un- 
dergo two-photon decay to l 1 ^ all have "1+1" decays in 
which the atom first undergoes an allowed single-photon 
emission to an intermediate level, and then undergoes a 
second allowed one-photon emission to reach the ground 
level: n l L -> n' 1 P° -> l 1 ^ (where 2 < n' < n). These 
1+1 decays are automatically included in the calculation 
of the two-photon spectrum using Fermi's Golden Rule 
[we will see this explicitly in Eq. lfl3|) ] and they turn out 
to dominate the net two-photon rate. In order to in- 
clude two-photon transitions from n > 3 levels in the 
multilevel atom code, we need to distinguish "true" two- 
photon decays from 1+1 decays. It is sometimes said 



that in a two-photon decay the two photons are emitted 
"simultaneously," but one must be careful in making this 
statement because the uncertainty principle dictates that 
one cannot measure the time of emission of the photons 
more accurately than the reciprocal of the frequency res- 
olution. Rather, one must return to the physical picture 
of recombination and remember that rare processes such 
as two-photon decay are potentially important because 
the He I resonance lines have a high optical depth and 
hence a high probability of re-absorption of any radia- 
tion emitted in those lines. In contrast, photons emitted 
outside of the resonance lines have a low probability of 
re-absorption (unless they later redshift into a line). 

Based on this picture, we can construct a "practical" 
definition for two-photon decays as follows: radiation 
emitted farther than some arbitrarily specified distance 
Ai/ cut from the nearest He I resonance line will be said 
to originate from a "nonresonant two-photon decay," and 
radiation emitted within Az/ cut of a resonance will be 
said to originate from a "resonant" or "1+1 decay." The 
nonresonant decays exhibit a continuous spectrum and 
can be treated in the same way as two-photon decays 
from 2 1 S. The full emission spectrum of the resonant 
decays is not identical to the usual Voigt profile, and 
the differences will have to be treated by modifying the 
line radiative transfer analysis. We will consider nonres- 
onant two-photon decays in Sec. |IV] and resonant decays 
in Sec. [Vj the corrections to the recombination history 
turn out to be small in both cases. This section will be 
concerned exclusively with obtaining the rate coefficients 
for two-photon decay, which we will find to be much less 
than estimated by DG05 across most of the two-photon 
spectrum. This is the reason why we find only a small 
correction from the n > 3 two-photon decays whereas 
DG05 found an effect of several percent in x e . 

We will also consider Raman scattering from the ex- 
cited levels to the ground level; the two processes, while 
physically distinct, are related by crossing symmetry and 
hence share many characteristics, including the existence 
of 1+1 resonances and the associated subtleties. 

The outline of this section is as follows. The formu- 
las for two-photon decay and Raman scattering in quan- 
tum electrodynamics are introduced and summarized in 
Sec. IIII Al The DG05 estimate for the rate coefficients 
is recalled in Sec. IIII B) and in Sec. IIII Cl we explain why 
their rates are too large for the high n levels. Finally, 
Sec. IIII Dl presents our calculation of the two-photon de- 
cay rates, which are much less than those of DG05, except 
near resonance. 



A. Rates 

On account of the electric dipole selection rules, two- 
photon transitions to the ground state of He I are al- 
lowed only from spin-singlet levels with even parity and 
L e {0, 1, 2}. In two-electron atoms only doubly excited 
levels can have L = 1 and even parity, and these are 
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inaccessible at recombination-era temperatures (they lie 
~ 60 eV above the ground state, whereas the ionization 
energy is 24.6 eV); thus we restrict our attention to the 
v}L levels with L £ {0, 2}. Also, the two photons emerge 
with frequencies v and v' that satisfy the energy conser- 
vation condition 

v , yl= E(n^L)-E(l'S) _ AEjn^L) (n) 

The two-photon decay rate from the n 1 L level to the l 1 ^ 
level of helium is then given by 



dr a 6 v 3 v' 3 



dv 108(2i + 1)11' 



6 [l+M(v)][l+M(v')}\M 2l \ 2 , 



where 7Z = 3.29 x 10 15 Hz is the Rydberg in frequency 
units, and the dimensionless amplitude is 

X 27 = ao 3 ^(l 1 5||d||n' 1 P )(7i' 1 P°| IdH^L) 



1 



^AEin^P ) - hv 

1 > 
AEin^P ) - hv' j 



(13) 



where d is the electric dipole moment operator and we 
have used cgs units. [Eq. (fl2|) is equivalent to Eq. (59.28) 
of Ref. [HI after appropriate manipulation of reduced 
matrix elements.] Note that the summation here is over 
continuum levels with 1 P° symmetry as well as discrete 
levels. The total two-photon decay rate is 



i r AE 
T 2l (n 1 L^l 1 S) = - J 



—dv, 
dv 



(14) 



where the factor of 1/2 occurs because we count each 
decay twice by integrating over the whole spectrum. 

One can see that the amplitude M. posesses a pole at 
each frequency v corresponding to an intermediate n' 1 P° 
level. Correspondingly, there is a branch cut (i.e. a con- 
tinuous distribution of poles) for frequencies correspond- 
ing to the n' P° continuum. Since the rate T oc \M\ 2 , the 
poles give rise to resonances in the cross section, which 
have the characteristic oc v~ 2 structure. As is usual 
in quantum mechanics (e.g. Sec. V§18 of Ref. [33|). 
the total rate is rendered finite by giving the energies 
E{n' 1 P°) a small imaginary part E — > E + iT/2, where 
r is the width of the state. The imaginary part changes 
the cross section in the resonance to the characteristic 
Lorentz form (which becomes a Voigt profile in the co- 
moving frame due to thermal motion of the atoms) . The 
resonances at < hv < AE(n l L) give rise to the al- 
lowed decays where the atom decays from n 1 L to n' 1 P° 
by emission of a single electric dipole photon, and then 
proceeds to decay to 1*S by emitting a second photon. 
These "1+1" decays are in fact not distinct physical pro- 
cesses from two-photon emission. Rather, the damping 



wings of the lines from 1+1 decays merge continuously 
into the two-photon continuum. 

A phenomenon related to two-photon decay is Ra- 
man scattering from n l L to l 1 ^ through an intermediate 
n' P° state. This has the same selection rules as two- 
photon decay. If the incoming photon frequency is v and 
the outgoing frequency is v' , we have 



AE{n r L) 



(15) 



and the scattering rate (in number of scatterings per 
atom in the n}L level per second) is 



(12) -T, 



dr a 6 v 3 v' 3 



dv 108(2L + 1)W 



-M{v)[l+N{v')]\M 



Raman | 



(16) 



Because of crossing symmetry, we may obtain the Ra- 
man scattering matrix element by analytic continuation 
of Eq. ifl3)l to negative frequencies: 



M Raman = « " E^ 5 ! I d l l"'^) ^ P °W A \ W L) 



XAEin^P ) + hv 

1 > 

' AE(n' l P°) - hv' j 



The total Raman scattering rate is 



Raman 



dv 



-dv. 



(17) 



(18) 



The summations in Eq. l|T3j) and Eq. ifTTj) are in gen- 
eral nontrivial as they depend on the helium wave func- 
tions. Accurate calculations are available only for the 
2 1 S level, which is the only singlet level for which two- 
photon transitions are the dominant mode of decay. In 
the case of cosmic recombination however, the "block- 
ing" of allowed one-photon electric dipole decays by high 
line optical depth means that subdominant decay modes 
of highly excited states can become significant, and es- 
timates of their rates are required. DG05 was the first 
paper to consider these two-photon decays in the con- 
text of the cosmic recombination, and they introduced a 
simple scaling argument for the rates. We revisit the is- 
sue here and conclude that the decay rate is significantly 
smaller. 



B. DG05 estimate 

This section reviews the derivation by DG05 of the 
two-photon rate from highly excited states in He I. We 
present the key points of the derivation in the notation 
of this paper in order to highlight the most important 
assumptions in their paper and how they differ from a 
more detailed treatment. 



DG05 noted that the dipole matrix elements of the 
form {n' 1 P°\ |d| ^L) are largest for n' = n. In particular, 
in the limit of hydrogenic wavefunctions they show that 
for large n, 



K^LWdw^p ) 



_9_ 

To 



^kn^HdlKV) . (19) 



Therefore DG05 argued that far from the 1+1 resonances 
in the two-photon rate, the matrix element M.2~< should 
be dominated by the n' = n term. They also noted 
the near-degeneracy of the n 1 P° and n 1 !/ levels. If one 
treates this degeneracy as exact, and keeps only the n' = 
n term in the sum, one can show that the vacuum decay 
rate is 




AE^P ) 



where the ' DG ' superscript indicates that the DG05 ap- 
proximation is being used. The frequency integral is a 
polynomial, and hence is trivially performed. It results 
in a total decay rate of 

AP*»oc(2L + l)-i ( AE{nlL) 

x K^F'lldllra^L)! 2 Kn^Udlll 1 ^! 2 . (21) 

Note that within the DG05 approximation all the two- 
photon spectra are scaled versions of each other. DG05 
thus used Eq. (|2~T]l to re-scale the H I 2s — > Is decay rate 
of 8.2 s _1 to the highly excited levels in hydrogen and 
helium, i.e. they re-scaled A in proportion to the squares 
of the dipole matrix elements and the fifth power of the 
energy difference. This leads to the result (using hydro- 
genic values for the (n 1 P°\ |d| |n 1 L) matrix elements) 



J2 (2£ + l)A (DG) (n x L) = 10540 



L=0,2 



n — 1 



n + 1 
lln 2 -41 



2/i 



(22) 



The most important result of this is the scaling for 
large values of n. At large values of n, AE{n}L) ap- 
proaches the ionization energy XHci, whereas the dipole 
matrix elements scale as (n 1 P°\ |d| \n}L) cx n 2 and 
(ri 1 P°||d||l 1 S') oc n~ 3 / 2 . This explains the large-n scal- 
ing of Eq. 4221): 



r£ G Vi)an. 



(23) 



This results in a very large contribution to the two- 
photon rate from large values of n. In fact, since the 



occupation probability of the large-n states approaches a 
constant as n — > oo, the total 2-photon decay rate to the 
ground state from highly excited helium atoms diverges 
as oc n oc n 2 in the DG05 approximation. DG05 cut off 
the sum at n ~ 40, since for larger n the "size" (~ a^n ) of 
the excited atom is comparable to the wavelength of the 
photon and hence the dipole emission formula is no longer 
valid. This nevertheless leads to a very large speed-up of 
He I recombination. 



C. Large n behavior 

Unfortunately, the simple approximation of taking 
only the n' = n term in the summation fails for large 
n. Indeed, it has been found for the highly excited states 
of hydrogen that the actual scaling of the two-photon de- 
cay rate is dk/dv oc ?i~ 3 [39]. Here we recall the physical 
argument why the scaling is n -3 , and then show that 
this arises due to a near-exact cancellation of matrix el- 
ements for large n. The argument has been given in a 
rather complicated and general form in Refs. |4Cl [UJ, 
however we present a simplified version here in order to 
highlight the key pieces of physics required, and see that 
the same argument applies to helium. 

Suppose we re-write the analogue of Eq. lfT3]) for H I 
in the form 

M 2l = a ~ 3 [(^)\\d\\nl) + <vf(z/)||d|H] , (24) 

where the states |\E' m (i/)) are defined by [12] 

|tf m (i>)) = [H- E{ls) - hv]- l dm\ls). (25) 

(Since the dipole operator d m is spin 1 it has three com- 
ponents m = —1,0, +1 and hence \^f m (u)) actually con- 
sists of three states.) Now the wave function of the | Is) 
state is localized near the nucleus, with an exponential 
falloff in the classically forbidden region. The dipole op- 
erator d m simply multiplies this wave function by a poly- 
nomial which does not affect the fact that there is an ex- 
ponential falloff. The state Y$> m (v)) is then determined 
by solution of the inhomogeneous Schrodinger equation, 



2m e 



-V 2 



E(ls) - hv 



*m(^;r) = d m ip ls (r). 



(26) 

This equation has been extensively studied in the context 
of the response of a hydrogen atom to electromagnetic 
radiation (e.g. Ref. [43|). At large r, the source falls off 
exponentially. Since E(ls) + hv < 0, the operator on the 
left-hand side of this equation takes the form of a wave 
equation with imaginary wave number ( u k 2 < 0") at large 
r. Therefore its solution is exponentially decaying with 
r. Now we know that as n — ► oo, the wave functions of 
\nl) states near the origin approach the solution of the 
Schrodinger equation at zero energy: 



(small r), (27) 
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Radial matrix element function, f cvc (fi) 




FIG. 2: The function /(fi) appearing in Eq. I|28|) . which de- 
scribes the radial matrix elements for large n. Note that /(fi) 
is largest for small £2, implying that {n 1 L\ |d| \n n P°) is largest 



for small n 



The points mark integer values of Q.. 



so that the ipniirf's near the origin are all solutions to 
the same linear homogeneous differential equation with 
regular boundary conditions at r = 0. Therefore they 
are scaled versions of each other, with the normalization 
determined by the condition J \ip n i\ 2 d 3 r = 1. The nor- 
malization integral is dominated by regions with r 3> ao 
and it is well-known that it enforces ip n i oc n~ 3 / 2 in the 
large n limit. Given that ip n i oc n -3 / 2 near the origin, 
and that * m (^; r) is only significantly different from zero 
near the origin, it follows that M.2-y oc n~ 3 / 2 and that 
dT^'y/dv oc n~ 3 . 

Essentially the same argument applies to helium. Like 
the | Is) state of hydrogen, the 1 1 1 S") state of helium is 
concentrated near the origin with an exponential falloff 
as either electron is moved to large distances from the 
nucleus. 

Mathematically, the only way to reconcile the .M27 oc 
n~ 3 / 2 scaling from the above argument with the oc n 1 ^ 2 
scaling obtained by including only the n' = n intermedi- 
ate state is that there must be a near-exact cancellation 
of contributions to .A/f 2 7 - We can show that this will hap- 
pen by examining the behavior of the matrix elements for 
large n and small s = nl — n. This is considered in Ap- 
pendix El where it is shown that 

(n^Udlln'V) a {-l) L > L 1 ^ 2 (-l) s ea n 2 f cyc (s + S 1L ) 

(28) 

(c.f. Eq. IA17|) . where 5\l is the differential quantum de- 
fect (<5io = 0.152 and Su = 0.014) and / cyc is the Fourier 
transform of the cycloid function (Eq. IA18|) . [The ex- 
istence of an asymptotic limit of the form in Eq. |28|) 
appears to have been noticed by Refs. [H, [45|, and the 
analytic expression for f rvr was derived, albeit in a dif- 
ferent form, by Ref. [46|.] The function / cyc is plotted 
in Fig. [2 What is of note here is that the matrix el- 
ements with s equal to a few are of the same order of 
magnitude as those with s = (n' = n). Therefore one 
should include them when obtaining the matrix element 
M.2-I- Since the matrix element (l 1 S , ||d||n' 1 P°) scales 



asymptotically as n'~ 3 ^ 2 , and the energy AE{n' x P°) ap- 
proaches a constant at large n' , these can be considered 
constant for \s\ <C n. We may thus include the values 

With S m \n < S < St, 



M 



2" 



"0 



by writing 
(l 1 S||d||n 1 .P )(-l) £ >.L 



1/2 

> 



E 



(-l)7cyc(* + <Jn) 



AS(ri 1 P ) - hv 



1 



AE^P ) - hv' 



(29) 



The first line in this equation scales as ti 1 / 2 , which when 
squared gives the DG05 scaling dh.2 1 /dv oc n. One must 
be mindful of the second line however, which modifies the 
prefactor of n 1 / 2 in the asymptotic scaling of .M27- One 
would expect to get a better estimate of the asymptotic 
scaling by taking the limits s mnl — ► — 00 and s max — > 00. 
However we show in Appendix [A] that (see Eq. IA19|) 



E (-ir/cycC*- 



0. 



(30) 



Therefore for large n one expects the contribution to .M27 
from states with n! near n to be n 1 / 2 times something 
approaching 0. This of course implies that one cannot 
find the large-n behavior of .M27 by the DG05 argument 
(except possibly by considering higher-order corrections 
to Eq. [28]) - one can only say that it scales slower than 
n 1 / 2 . The near-cancellation is illustrated graphically in 
Fig. [3] for the 25 ^ -> l 1 ^ decay. 

In summary, we conclude that (i) the actual large-n 
behavior of the two-photon decay rate is dA/dv oc ?i~ 3 , 
and (ii) the apparent discrepancy between this and DG05 
is due to a cancellation of the matrix elements as sum- 
marized by Eq. ([50]). 



D. Rate estimates 

We have estimated the two-photon transition rates for 
small n by direct summation of the matrix element prod- 
ucts in Eq. (fl3|) . We note that such a calculation does 
not require a detailed re- analysis of the atomic physics, 
at least to a first approximation, because all one needs 
to know are the energies and dipole matrix elements for 
the 1 S, 1 P°, and 1 D levels, which have already been cal- 
culated. (The exceptions are the continuum levels, for 
which the information available is more sparse, neverthe- 
less as we argue below there are detailed calculations for 
the first several eV of the continuum, which are domi- 
nant.) We have obtained these as follows: 

1. For the n 1 S^n' P° with n,n' < 9, we use the os- 
cillator strengths from Ref. [48l |. 
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Helium matrix elements from n=25 
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10" 2 <25lS||d||n';P°>/ea — 
10 <1 1 S||d||n ,1 P°>/eaQ * 
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FIG. 3: The dipole matrix elements |d| |n' 1 J 50 ) and 

(n 1 5'JJd||n' 1 P°) for n = 25, computed using the formula of 
Ref. [471 ] (for l 1 ^) and the Coulomb approximation (for 25 1 S). 
Note that the matrix elements with n' close to, but not equal 
to, n are comparable in magnitude to {r^S] |d| |n 1 P°) but have 
negative instead of positive sign. Therefore these contribu- 
tions can interfere destructively with the n 1 P° intermediate 
level in the matrix element .M27. 



Helium n=2 2-photon decay rates 
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FIG. 4: The 2-photon spectrum from the 2 1 S level of He I. 
Note that there are no resonances in the spectrum since there 
are no energetically allowed 1 + 1 electric dipole decays from 
2 X S. The points are the more detailed calculations by Drake 



2. The l 1 ^ - 10^° oscillator strength is from 
Ref. Hi]. For ^S-n^P" transitions with n' > 10, 
we have used the asymptotic formula for the oscil- 
lator strength from Ref. [47j |. 

3. For the n 1 S-n' 1 P° transitions with n,n' > 10, we 
used the Coulomb approximation [HI]. 

4. For transitions from r^S to the 1 P° continuum, 
we used the TOPBase photoionization cross sec- 
tions (50I ] converted to matrix elements in accor- 
dance with 



4a dn' 
^ hV dE 



[n l S\\d\\n' X P°) 



(31) 



This turns the continuum contribution to Eq. fl3j) 
into an integral over energy. 

5. Dipole matrix elements for S-P° transitions are ob- 
tained from the standard formula, 



[n 1 SW&Wn' 1 P°) 



3e 2 alUf n i S ^ n apo 
Ein^P ) - E(n 1 S)' 



(32) 



where dn'/dE is the density of continuum states. 



The helium atom wavefunctions for m = are all 
real and hence (n 1 »S , ||d||r7,' 1 P ) is purely real, how- 
ever a sign ambiguity exists. We have taken the 
sign for S-P° to be negative for n 7^ n' and posi- 
tive for n — n', as this is what is found using the 
Coulomb approximation or hydrogenic wavefunc- 
tions. 

6. The P°-D dipole matrix elements are taken to be 
hydrogenic. 

Note that this approach is expected to break down for 
continuum 1 P° levels with very large energies. In partic- 
ular the continuum l P° wave functions become less hy- 
drogenic at higher energies where the outer electron pen- 
etrates deeper into the He + "core," and it becomes very 
non-hydrogenic as one approaches the double-excitation 
resonance region 60 eV above the ground state. At still 
higher energies there are multiple continua, so it is no 
longer valid to compute bound-free matrix elements us- 
ing Eq. (|3"Tj) - the matrix elements actually contain in- 
formation that is not contained in the cross section. For- 
tunately, these subtleties have little effect at the level of 
accuracy required here: we find that neglecting contin- 
uum levels with energies more than 0.57?. (6.8 eV) above 
threshold makes at most a change of 30% (3 X 5 — > l x 5) 
or 1% (3 1 D — ► l 1 ^) to dkjdv, except in the immediate 
vicinity of the nulls (whose positions are slightly shifted) . 
Since we will find that the total correction to the recom- 
bination history due to nonresonant two-photon decays 
is ~ 4 x 10 -4 , we believe that our basic conclusion that 
nonresonant two-photon decays are unimportant is ro- 
bust even if the rate estimates are off by several tens of 
percents. 

We show the two-photon decay rates we obtain for the 
He I 2 1 S level in Fig. |4j The total decay rate we obtain 
is 49 s -1 , in comparison with the more detailed atomic 
physics calculations, which give 51 s _1 [13] ■ This provides 
a check on the accuracy of our method. 

Also of interest are the two-photon rates from the 
n = 3 (Fig. [5]) and n — 4 (Fig. [6]) levels. These two- 
photon spectra show resonances at the positions of al- 
lowed 1+1 transitions. We have shown the results using 
only the n' = n terms as a series of points in each plot; 
one can see that this is a poor approximation across most 
of the spectrum. In particular, far from the resonances, 
this substantially overestimates the rate because it ne- 
glects destructive interference between different interme- 
diate levels. 
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Helium n=3 2-photon decay rates 



1000 



'w 





3, 1 S-lls 




- 


3 1 D-1 1 S 






\ 




f 




31 



5 10 15 20 

Photon energy E, eV 

FIG. 5: The 2-photon spectrum from the 3 1 S and 3 1 D lev- 
els of He i. Note the resonance at 21.2 eV corresponding to 
the 2 1 P° intermediate level. There are also resonances at 
much lower energies corresponding to the optical transitions 
in He I, 3 1 S~2 1 P° and 3 1 D-2 1 P° . Also note the nulls in the 
two-photon rate from 3 1 5'. The series of points are the re- 
sults considering only the n = n term in the matrix element, 
Eq. (I13II . for 3 1 S (upper series) and 3 1 D (lower series). Note 
that keeping only this term is a poor approximation except 
at the very ends of the spectrum. 



Helium n=4 2-photon decay rates 
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FIG. 6: The 2-photon spectrum from the A 1 S and A 1 D levels 
of He i. There are two pairs of resonances corresponding to 
the 1+1 decays via the 2, 3 1 P° intermediate levels. The series 
of points are the results considering only the n' = n term in 
the matrix element, Eq. lfl~3"|) . for 4 S (upper series) and 4 D 
(lower series). 



IV. EFFECT OF NONRESONANT 
TWO-PHOTON TRANSITIONS 

Now that we have obtained the two-photon rates, we 
would like to understand how much He I recombination 
is modified by including them. The main contribution 
comes from the lower values of n, both because of their 
faster rates and because the lower-n states have higher 
occupation probabilities. This section considers the ad- 
dition of nonresonant two-photon decays from 3 < n < 5 
and nonresonant Raman scattering from 2 < n < 5, and 
finds a negligible effect. 



There is one subtlety involved in including higher-order 
two-photon transitions, which was recognized already in 
DG05. It is the existence of the 1 + 1 resonances, which 
cause the two-photon rate to be very large when the pho- 
tons are emitted in allowed electric dipole lines. Photons 
emitted in these lines (i) have a high probability of being 
re-absorbed, and (ii) are in any case already included in 
the treatment of Paper I, which included all of the one- 
photon transitions. In this paper, we will handle this 
issue by dividing the two-photon spectrum into nonres- 
onant and resonant pieces, which are treated separately. 
Here "nonresonant" simply means that the emitted pho- 
tons are detuned from the 1+1 resonance by some mini- 
mum frequency offset Af cut . The idea is to show in this 
section that the nonresonant transitions have no signifi- 
cant effect on He I recombination, and then in the next 
section consider whether the approximations made in Pa- 
per I about resonant two-photon transitions are valid. 
The offset A^ cut is arbitrary and was chosen so as to 
make both the arguments in this section and the follow- 
ing section valid. Precisely the same subtlety arises in 
considering Raman scattering, which has resonances such 
as 2 1 S — > 2 1 P° — > l x S, and we handle the problem in 
precisely the same way. The choice of the frequency off- 
set that we use is O.Un^TZ = 460?i~ 3 THz for the offset 
from the He I n 1 P°-l 1 S line; the motivation is that we do 
not want our definition of "resonant" photons to overlap 
with the intercombination line He i] n 3 P°-l 1 S. (There 
is an overlap with the quadrupole lines [He i] r^D-^S, 
however as we argue in Paper III, these lines do not mat- 
ter anyway.) 

In the absence of the spectral distortion, nonresonant 
two-photon transitions and Raman scatterings can be 
trivially included in a level code as an additional rate, 

x (x n , L -^x lls e- E n^A ,(33) 
V 9i l s J 

where the term with xii$ accounts for thermal re- 
excitations of ground-state helium atoms determined via 
the principle of detailed balance. The two-photon and 
Raman scattering rates, T2j{n 1 L) and rR aman (n, 1 L), are 
obtained by integration of Eqs. lfl4|) and (fl8|) with black- 
body radiation profiles, except that regions in the integral 
where the higher-energy photon lies within A^ cut of an 
allowed resonance are excluded. We have obtained fit- 
ting formulas for the two-photon rates, which are given 
in Table |H We have also included nonresonant Raman 
scattering from the 2 X S level, which is well fitted by 

r Rama „(2 1 ^ = U.St^e- 2 - 125 /' a" 1 , (34) 

where t = T r /4000K is in the range 1 < t < 2. For 
1 < t < 2, Eq. (|34j) and the formulas in Table U agree 
with our numerical calculations to within 1%, which is 
probably better than the accuracy of our rates. 

Fig. [7| shows the change in the electron abundance due 
to nonresonant two-photon transitions. The effect is at 
the level of a few times 10~ 4 and can be neglected. 
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TABLE I: Fitting formulas for the two-photon rates 
r 2 ~ / (n 1 L) +r Raman (n 1 L) appearing in Eq. ((33j) for the n = 3, 
4, and 5 levels of helium. The rate in units of s _1 is written 
as a + bt, where t = T r /4000 K. The formulas are valid in the 
range 1 < t < 2 of interest to helium recombination. Note 
that these numbers include only the nonresonant contribu- 
tion, defined as having the emitted photon at least 0.14n' 3 TZ 
from the He I n' 1 P°-l 1 S line. 



Upper level 


a 


b 


Upper level 


a 


b 


3 1 S 


20.5 


16.0 




94.0 


8.0 


4 X 5 


10.2 


12.5 


A 1 D 


42.3 


6.0 


5 X 5 


6.1 


9.2 


5'D 


21.5 


5.0 



Ax e from nonresonant 2y decays from n>2 levels 




-5x1 0" 4 1 1 1 1 1 

1600 2000 2400 2800 3200 

z 

PIG. 7: The contribution of the non-resonant part of the two 
photon rate in He I for n > 2 produces a maximum change in 
the free electron fraction of severalx 10~ 4 . 



V. THE 1 + 1 RESONANCES AND FINITE 
LINEWIDTH 

In Sec. IIV| we considered the influence of the non- 
resonant two-photon transitions on He I recombination. 
We know however that the total two-photon transition 
rate is dominated by the 1 + 1 resonances (except in the 
case of the 2 1 S level, which has no such resonances). If 
this additional rate is naively added to the recombination 
equations in the manner of Eq. (|33| , He I recombination 
becomes described by the Saha equation. However we 
know that this naive addition is incorrect because pho- 
tons emitted within resonance lines with lower level l 1 ^ 
will likely be re-absorbed. In order to understand the ef- 
fect of resonant two-photon transitions, we must under- 
stand the transport of radiation within the He I n 1 P°- 
1 X S lines. We presented a simplified analysis of this in 
Paper I, where photons were injected into the line by res- 
onant two-photon emission and H I recombination, trans- 
ported by coherent (Rayleigh) scattering and Hubble red- 
shifting, and finally removed by resonant two-photon ab- 
sorption and H I photoionization. The analysis in Paper 
I makes the approximation that the He I line is infinites- 



imally thin relative to variation in the radiation phase 
space density and phase space factors. The purpose of 
this section is to test the validity of these assumptions 
in certain special cases and understand the errors intro- 
duced. The basic method here is to reconsider the 2 1 P°- 
l 1 ^ line including the deviation from Voigt profile in the 
far damping wings, and including the deviation of the 
radiation profile from steady state. We incorporate these 
corrections into the level code and show that the modifi- 
cation to recombination is small (|Ax e | ~3x 10 -4 ). 

The specific assumptions made in Paper I that we 
would like to test arc: 

1. The two-photon emission profile can be described 
by a Voigt distribution, i.e. we neglected the 
possible interference with neighboring 1+1 reso- 
nances, and the variation of the photon phase space 
factor v z v' 5 and the phase space density factor 
[1 +M{v))[l +N{v')] (cf. Eq. 021) across the line 
width. (A similar assumption applies to our treat- 
ment of Raman scattering and two-photon absorp- 
tion.) 

2. The He I line was treated as being in steady state, 
i.e. we assumed that the rate of injection of photons 
equaled the loss rate. In reality, there are always 
a few photons within the line, and as this num- 
ber of photons increases (or decreases) there is a 
corresponding speed-up (or slow-down) of He I re- 
combination. 

We will examine these assumptions here in the con- 
text of the 2 1 P°-\ 1 S line, which was found in Paper I to 
produce the most important effect. To simplify the calcu- 
lation, we will also assume when calculating line shapes 
that the excited levels in He I are in equilibrium. (This 
was found in Paper I to be a good approximation and 
is described quantitatively in Paper III.) We will intro- 
duce the notation v- and v+ to denote the minimum and 
maximum frequencies of the resonance, i.e. 

v± = ^ii s _„ipo ± O.Un^TZ, (35) 

where here n = 2. 



A. Finite linewidth 

The calculations involving transport and incoherent 
scattering in Paper I made use of an approximate sym- 
metry of the thermal radiation field in the neighborhood 
of the line. Here, so long as the linewidth is negligible 
compared to k^T^/h, differences in the thermal radiation 
field on either side of the line can be neglected. In this 
section we will elaborate on this and argue that: 1) a 
linewidth of much less than kftT T /h means that a pho- 
tons are just as likely to be absorbed on either the red or 
the blue side of the line, and 2) that introducing a finite 
width to the line means that photons will be scattered 
differentially depending on whether they are on the red 



11 



1 1 S 2'P° 3 1 D 2 1 P° 1 1 S 




15B4 1.6678 1,6078 1584 



FIG. 8: An example of an incoherent scattering process in 
He i considered in Sec. IV Al 

or the blue side of the line. This will result in photons be- 
ing "pumped" redward by the asymmetry in the thermal 
radiation field across the line. 

Incoherent scattering is inherently a multi-photon pro- 
cess. Here, a photon is scattered and another is re- 
emitted and distributed over the line's profile, with no 
memory of the incoming energy. The only way for this 
reaction to proceed is for some number of other parti- 
cles to recoup the change in energy (as compared to co- 
herent scattering, where the photon's energy is exactly 
conserved in the atom's rest frame). 

Consider the case of incoherent scattering off He I l 1 ^ 
through 2 1 P° with an excursion to 3 1 D. Here the incom- 
ing photon (A) in the He I A584 line excites the atom to 
the 2 1 P° level. This excited atom resonantly scatters a 
second photon (B— >C) in the He I A6678 line via the 3 1 D 
resonance. Finally the atom decays back to the ground 
state, emitting a photon in the A584 line (D). This pro- 
cess can be viewed as a resonant two-photon absorption 
of photons A and B, followed by two-photon decay emit- 
ting C and D. 

In principle photon B could be any photon drawn from 
the blackbody radiation field, however because of the nar- 
row 3 1 D resonance in helium, the photon absorbed will 
almost always have energy Eb = AE(3 1 D) — Ea- The 
excited He I 3 X D atom then undergoes two-photon decay 
to the ground state via the 2 1 P° intermediate level (i.e. 
it emits photons C and D). To a very good approxima- 
tion, the energy distribution of D is independent of Ea 
- hence the term "complete redistribution" - and in the 
vicinity of resonance it has the form of a Lorentz pro- 
file (or a Voigt profile in the comoving instead of atom 
frame) . 

Now, suppose that the atom absorbs photon A on the 
blue side of the A584 line. Then it can be absorbed in 
combination with a photon B on the red side of A6678, 
whereas if photon A is on the red side of the A584 line 
then it requires B to be on the blue side of A6678. Since 
in a blackbody distribution for photon B there are more 
photons on the red side of the line, this means that there 
is an enhancement in the cross section for absorbing pho- 
ton A from the blue side of A584, and a suppression for 
absorbing it from the red side of A584. This means that 
(even in the absence of Hubble redshifting) A584 photons 
spend on average more time on the red side of the line, 
so that Af is greater there. 

The same conclusion could also have been reached by 
a thermodynamic argument: since incoherent scattering 




FIG. 9: Representation of complete redistribution as a two 
photon process, with one photon from the thermal distribu- 
tion. In the left frame, a photon is absorbed on the red side 
of the line, then assisted to 3 X D by a black-body photon. In 
the right frame, a photon is absorbed on the blue side of the 
line and assisted by a lower energy black-body photon. The 
virtual levels have energies offset from E(2 P°); the fractional 
difference between the forward and backward scattering rates 
is of order the frequency difference times h/ftsT?. Because 
there are more low-energy thermal photons, scattering to the 
blue side of the line is slightly preferred. 



changes the energy of the A584 line photons by exchang- 
ing their energy with that of the A6678 photons (and with 
other low-energy photons if we consider the other lines 
connecting He I 2 1 P° to other excited levels), and the 
radiation in these lines is essentially blackbody, it follows 
that photons near the A584 line will then be driven to- 
ward a Bose-Einstein distribution with temperature T r 
and some chemical potential determined by the total 
number of such photons. Since Af <C 1, this is equivalent 
to a Boltzmann distribution, N oc e -hAv/k B T t ^ \y e will 
see this behavior mathematically from Eq. (|39|) . (Note 
that whether Af oc e ~ hAv ' kBT ' is actually achieved de- 
pends on whether incoherent scattering can operate effi- 
ciently before Hubble redshifting moves the photons out 
of the line, a question that can only be settled by solving 
the equations.) 

In summary, we have argued that incoherent scatter- 
ing through a finite linewidth enhances the phase space 
density on the red side of the A584 line. In the limit that 
the width of a line is taken to be negligible compared to 
knTt/h, incoherent scattering redistributes photons and 
pushes the radiation phase space density near line center 
to some constant A/l, in equilibrium with the line. This 
flattening tendency is implicit in the analysis of Paper I, 
and has been noted several times in the recombination 
literature HJUl]. 

Viewed in this way, incoherent scattering is the sum 
of two-photon scattering processes for which the excited 
level is an intermediate state (resonance) in the full two- 
photon rate. The goal, then, is to consider the full ex- 
pression for all two photon processes, and separate the 
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transport physics around one of the intermediate state 
resonances, e.g. 2 1 P°. It is possible, then, to write an 
effective one-photon transport equation (where the other 
photon is drawn from the black body) for incoherent scat- 
tering to this intermediate state. Aside from the change 
in the line profile, the finite linewidth introduces two new 
pieces of physics: the tendency to drive the radiation 
spectrum to AT oc e - hv / k BT r ms t e ad of a constant, and 
the fact that the line is not exactly in steady state (i.e. 
there is a dAf/dt term in the transport equation). This 
section considers these issues and introduces a crude cor- 
rection to the rate equations. We incorporate this in the 
level code and find only a small correction (a few times 
1CP 4 ). This correction is not included in the final version 
of the recombination history presented in Paper III. 



Line transport with complete redistribution and 
no H I opacity 



The case we consider here is that where the H I 
opacity within the line and frequency diffusion due to 
Doppler shift in repeated resonant scatterings can be ne- 
glected. This is useful for testing assumption ^43] on our 
list (Sec. [Vj. The assumption of negligible H I opacity 
is valid in the early stages of He I recombination, i.e. 
2200 < z < 2800. The frequency diffusion was included 
in Paper I and neglecting it was found to introduce no 
significant error: \Ax e \ < 2 x 10~ 4 . 

The resonances in consideration are optically thick 
and the radiation rapidly approaches equilibrium around 
their line centers. Because of this, Doppler broaden- 
ing can be neglected to a good approximation and we 
can consider the radiative processes as occurring in the 
atom's rest frame. The differential equation describing 
the radiation field is 



dJV „ dAf . c 3 n H ^ dT t 



uj\ uj\ v /in v - ^ "a i 

~dt = + W l^ Xi ~^J 

i 

E(27+Raman) / -. , r 
<x\ {v)N, 



(36) 



where the sum is over excited levels of He I that can 
undergo two-photon decay or Raman scattering to the 
ground state, dYijdv is their rate of producing line pho- 
tons per unit frequency, and -( 2 T+ Raman ) j s th e cross sec- 
tion for removing line photons via two-photon absorption 
or Raman scattering to level i. This is a strong function 
of v, but we will drop the explicit argument to stay con- 
cise. (Though officially a three-body process, it is pos- 
sible to define a cross section for two-photon absorption 
of a 2 1 P°-1 1 S line photon since the other photon comes 
from much lower energies where the CMB can be treated 
as a blackbody.) By detailed balancing of the level i 
contributions to the second and third terms on the right 



hand side, we find 

C 3 n H dTj gj ^-AEJkyT, _ n co .(27+Raman) e _/ u ,// !B T t 
dv gil g 

(37) 

which allows us to derive the cross sections for two- 
photon absorption to each level. Then since the excited 
levels are in equilibrium with each other we have 



9i 



92 1 P° 



-x 2 ip°e 



-[Ei-E(2 1 P°)]/k s T T 



(38) 



Combining Eqs. 1(55)1 ; (|37)l . and {381), we S et 



dt dv 

— TT-HCXiisC 



(27+Raman) 



JSf _ X21p ° c -hA V /k B T, 
3a;iis 



(39) 



There are several approaches available for solving 
Eq. f39)) . We will take an approach that allows us to 
separate the effects of the line profile from the steady 
state approximation. The method is to multiply the left 
hand side of Eq. ((39|) by an artificial expansion parame- 
ter e, which will eventually be taken to equal 1. We may 
then expand 



Af = Af + eA/k + e 2 A/2 + 



(40) 



equating coefficients of in Eq. ((39]) then leads to the 
following situation. For j — 0, we find 



= His 



dAfp 
dv 



-7lHCXiis<J 



(27+Raman) 



x (Af -Af h e- hAv / k * T *), 



(41) 



where A/l = x<pp° /Zxiig. That is, Afo satisfies the 
steady-state equation. The higher-order terms satisfy 



dM 



_jti = Hv^L n^sa^—W, (42) 
at ov 

for j > 1. Since photons enter from the blue side of 
the line, the boundary condition Af(v+) is satisfied; the 
Taylor expansion of this condition in e is that Ao(^+) = 
Af(v+), and Afj{v + ) = for j > 1. We may think of 
the Afj for j > 1 as successive corrections to the steady- 
state solution. For each j, a numerical solution may be 
obtained by starting at v = v+ and using a stiff ODE 
integrator in the redward direction until we reach U- . 

In order to translate our results for the line profile into 
effects on recombination, we need two numbers. One of 
these is the photon phase space density Af{vJ) emerging 
from the red side of the line, necessary to compute feed- 
back. The other is the net decay rate to the ground state, 
which is obtained by subtracting the downward from the 
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upward rates: 



EdTi 



8ttv 2 



(2*7+ Ram an) 



dv. (43) 



The downard and upward rates nearly cancel, so numer- 
ically the best way to compute this is not to evaluate 
Eq. (|43| directly from the solution, but rather to use 
Eq. lj36j) to re-write it as 



n hc 3 



(44) 



The steady-state solution is obtained in Eq. (|44]) by set- 
ting M = Afo and 6 = (i.e. dropping the time derivative 
term). The first-order solution in e is 



x i 



8tt^ 



Hv 



d(N + M) 3Au 



o>/ 



dv. (45) 



The line profile is shown in Fig. [10] for a typical set 
of parameters, and is compared with the infinitesimal 
linewidth approximation, the steady state solution, and 
the analytic model of Appendix [B] The most important 
property of the solution, which is generic, is that A/" > A/l 
for Av < 0. That is, the effect of using the full Ao +M 
in Eq. (|45]) instead of just a step function at the line is 
to enhance the decay rate and accelerate recombination. 
On the other hand, dNa/dt < 0, so the correction due 
to the line not being exactly in steady state is of the 
opposite sign: it delays recombination. 



C. Inclusion in the level code and the effect on 
recombination 

The basic strategy in including the finite linewidth ef- 
fects in the level code is to determine the corrections to 
the phase space density N{y-) on the red side of the 
line and the net downward transition rate x^. This sec- 
tion describes how we do this, and the results when the 
correction is incorporated in the level code. 

The recombination level code depends on 
and the reaction rates implied by finite 
linewidth. In general these depend on the parame- 
ters {z,XHei,xnei,Nh,Nh,N'(i / +),Af(v+)}. Since the 
equation for M is linear, A/"o depends linearly on the 
parameters {A/l, A/"(v + )}, and A/l depends linearly on 
the parameters {Al, A/l, A/"(v+), J\(v+)}. Also ±h g i 
enters only via Ao, which is the source for A/l (c.f. 
Eq. [42]). From this one can see that the phase space 
density may be written as 

N{v) = c (v)ATl + Ci(i/)A/l 

+c 2 {v)N{v + ) + c 3 (v)Af(i> + ) 

+:r H ci [c 4 (^)AA L + c b {v)N(v+)] , (46) 
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FIG. 10: The phase space density in the 2 1 P°-1 1 S line 
using the parameters z = 2175, xnei = 0.02825, A/L = 
6.307 x 10" 16 , M(y+) = 5.277 x 10~ 19 , i Ho i = 0.1596H, 
A/L = -2.444 x lQ- lh H, and M{v+) = -4.073 x lO" 17 ^. 
These parameters occurred during the first feedback iteration 
of the recombination history. The solid line shows the steady- 
state solution A/o, while the long-dashed line is the first-order 
correction A/o +A/i. The short-dashed line is the analytic ap- 
proximation to the steady-state solution from Eq. (IB 1 T|) : note 
that it is plotted only for Au < 0. 



where the Ciiy) depend on z, SHci, and cosmological pa- 
rameters. Thus if we want A/"(v_), then for each cos- 
mology an interpolation grid can be constructed to give 
Ci{v-) in terms of the independent variables z and iHei- 
A similar result holds for ±j since it is a linear function 
of Wo, A/o, and A/i- 

The easiest way to incorporate the new effect in the 
level code is actually to calculate the correction to N(v- ) 
and In the case of infinitesimal linewidth, no con- 
tinuum opacity, and high optical depth (literally, neg- 
ligible probability of a photon redshifting through the 
line without undergoing an incoherent scattering - see 
Appendix D of Paper I), we have P csc = t^ 1 . In this 
case, the photon phase space density on the red side 
of the line is A/l and the downward transition rate is 
87r_ff^ 1 3 inc (A/'L — A/+)/riHC 3 . If we ask about the photon 
phase space density at v_ < Un n e, and specify the incom- 



ing (blue-side) phase space density at v + 
becomes 



this 



Af(v-,Z-) — > A/l {z) and 
x l( z ) 



" c [A/l(z)-A^ + ,z + )], (47) 



where 



71HC' j 



l + z ± = J^-(l + z ) 

^linc 



(48) 
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since it takes a finite amount of time for photons to red- 
shift through the line. One may thus define a "finite 
linewidth correction" 

8N{v- ,z~)= Af(v- ,z~) — A/l (z) (49) 

for the phase space density on the red side of the line 
(used for feedback), and a similar correction 

Sx^z) = x^z) - ^ Hv f c W^z)-N{ V+1 z + )] (50) 

7T.HC 

for the transition rate. 

We have re-run the level code with Eqs. lj49j) and fBTjj) 
incorporated and "turned on" from z = 1500 to 3400. The 
change in x e is shown in Fig. [TTJ The correction is be- 
lieved to be most accurate for z > 2200 when continuum 
opacity is negligible. At lower redshifts, the corrections 
of Eqs. (@9j and ^ are not reliable. For Eq. (J49j) this is 
not a major deficiency because at these redshifts feedback 
[the only process affected by Af(i/-, z_)J is unimportant. 
For Eq. (|50)) there is an error introduced, however we 
expect that the change in x e at z < 2200 (when contin- 
uum opacity is significant) is small because it is only in 
the far damping wings that the corrections described in 
this section are significant, and continuum opacity makes 
the line center more important relative to the damping 
wings. (This is because continuum opacity allows pho- 
tons to be removed from the line center, whereas without 
continuum opacity photons can only escape the line by 
redshifting out of the red damping wing.) 

The modification to the recombination history result- 
ing from these changes is shown in Fig. [TTJ We see that 
the total effect reaches a maximum of 0.03% in the free 
electron fraction. This is much smaller than the other 
effects and comparable to other errors in the code, so we 
have made no attempt to correct for the deviation from 
Voigt profile or change in e -W fc BT r across the j me m 
the rest of this series of papers. Since the correction is 
of order the numerical accuracy of the code and involved 
such a major change to the treatment of the all- important 
2 1 P°-1 1 S resonance line, we do not claim that the de- 
tails of Fig. [TTJ are robust; rather we view the results only 
as confirmation that the effects considered are small. 



VI. DISCUSSION 

This was the second paper in a series devoted to cos- 
mological helium recombination. Here, we examined 
the problem of two-photon decays in He I, extending 
the standard treatment which only accounts for the de- 
cay from the 2 1 S level and ignores the effect of stimu- 
lated transitions and absorption of the spectral distor- 
tion. We also considered Raman scattering from excited 
levels in He I to the ground level (l 1 ^), an effect that is 
distinct from, but closely related to, two-photon decay. 
All of these effects change the electron abundance x e at 
the level of several hundredths of a percent at redshifts 
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FIG. 11: The change in the recombination history from the 
modified treatment of the 2 1 P°-1 1 5' resonance. Note that the 
effect on the electron abundance is very small: a few parts in 
10 4 . This fi gure should only be interpreted as an estimate of 
the magnitude of the correction (see text). 



z > 1800. This results in a change of similar magnitude 
in the C/s (the precise relation will be quantified in more 
detail in Paper III), which is negligible for cosmic He I 
recombination studies. 

Our findings regarding the significance of two-photon 
decays from the n > 3 levels of He I differ from some re- 
cent statements in the literature, most notably Dubrovich 
& Grachev who found a much larger effect. The 
main reason for the difference is that we find smaller two- 
photon rates dA 2l /dv because of destructive interference 
among different intermediate states in the two-photon 
amplitude (Eq. 1X3]) in most parts of the two-photon con- 
tinuum. An exception occurs in cases where the two pho- 
tons emitted in a two-photon decay are near an allowed 
"1+1" sequence of decays such as i^S -> 2 l P° -> l x 5. 
These 1+1 decays correspond to resonances in the two- 
photon decay rate at the frequencies corresponding to the 
one-photon lines (in our example, the 3 1 S — > 2 1 P° and 
2 1 P° — > l x S lines). This results in the total (frequency- 
integrated) rate being very large. This does not lead 
to a rapid speed-up of recombination however, because 
the photons emitted in the optically thick resonance lines 
have a very high probability of re-absorption. In order to 
complete the recombination calculation it is necessary to 
split the photons into resonant and nonresonant regions. 
The nonresonant regions are handled in the usual way for 
two-photon decays, i.e. they lead to an additional rate 
that is included in the rate equations. The two-photon 
decays in the resonant regions are treated as sequences of 
one-photon decays, with the two-photon effects leading 
to a modified line profile since with multiple intermediate 
states the Lorentz curve no longer accurately describes 
the line profile (in the atom's rest frame). It is essential 
in this analysis that the treatment of the resonant re- 
gion takes into account the fact that the line is optically 
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thick, otherwise unrealistically fast recombination would 
be obtained. 

The analysis presented in this paper was aimed pri- 
marily at helium recombination, however most of the un- 
derlying physics is the same for hydrogen recombination. 
There are two-photon decays from the n > 3 levels in H I, 
and their rates dA 2l /dv scale as n -3 [3§| for the same 
reasons described here. These rates also possess reso- 
nances at the frequencies corresponding to 1+1 decays 
such as 3s — > 2p — > Is. In general hydrogen recombi- 
nation matters more for the CMB power spectrum than 
helium recombination, and in particular Wong & Scott 
[29| have found changes in the C/s of several tenths of a 
percent using rates much smaller than those of DG05. A 
full calculation for hydrogen would use the two-photon 
spectra dA2~ ( /dv, which could be computed by the same 
methods used here, and take into account the modifica- 
tion of the Lyman line profiles due to two-photon cor- 
rections. Such a calculation is beyond the scope of this 
series of papers, but should be a high priority for the 
CMB community. 



the combination of the inner electron and nucleus (of 
charge +e). Except at small r, its radial wave func- 
tion R(r) thus satisfies the Schrodinger equation R"(r) = 
—k 2 (r)R(r), where 



k 2 (r) 



2m e 



E 



e 
r 



h 2 L(L + l) 
2m e r 2 



(Al) 



For E < 0, this equation possesses a classically forbidden 
region r > r max , where r max is the solution to fc 2 (r max ) = 
0. In the classically allowed region, the WKB solution for 
R{r) is 



R(r) = (-1) 



i-L-l 



N 



■. cos (p(r), 



(A2) 



where the (— l)" _i_1 factor is chosen by convention to 
make the wave function positive near the origin for A^ > 
(it has n — L — 1 radial nodes) and the radial phase is 



k(r) dr. 



(A3) 
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APPENDIX A: DIPOLE MATRIX ELEMENTS 
FOR LARGE n 

This Appendix evaluates the dipole matrix elements of 
the form {n 1 L\ |d| {n^P ) for large n and small s = n' —n 
using the Wentzel-Kramers-Brillouin (WKB) method. 
This approach is useful since the dipole matrix elements 
of this form are dominated by large radii where the WKB 
method works (it breaks down at radii of order ao or less) . 
Our goal is to demonstrate the near-exact cancellation 
of contributions to M 27 in Eq. lfl~3|) that we mentioned 
in Sec. IIIICI We note that WKB-type solutions to the 
Coulomb approximation wave function have been previ- 
ously used for several other applications [46L [5l| . The 
formula presented here is actually equivalent to the spe- 
cial case of Ref. in which the eccentricity of the or- 
bit goes to 1, however we provide a simplified derivation 
here in order to show the fastest route to the key result 
(Eq.EH. 

For large n, the helium atom can be treated by the 
Coulomb approximation in which the outer electron (of 
charge — e) moves in the Coulomb potential defined by 



The normalization constant is taken to be positive, and 
to enforce the condition J \R(r)\ 2 dr = 1. For small L, 
we have r max w e 2 /(—E), the classically allowed region 
extends down to r <C r max , and then (for small L) 



k(r) = ^■y / 2m e (r 1 - r max ). 
From this we find 

\R{r)\ 2 dr 



(A4) 



2 , A^ 2 rr " 



3/2 

max N 2 , (A5) 



k(r) 4V2m e e 

where we have replaced cos 2 <p(r) with 1/2 since we in- 
tegrate over many oscillations of the wave function. This 
gives 



R(r) 



(-1) 



(^r max )V2( rmax / r _ i)V4 



cosip(r). (A6) 



In order to compute radial matrix elements with these 
wave functions for small s, we need to consider the effect 
on the wave function of small changes in k 2 (r) resulting 
from changes in E and L. In general there will be a very 
small change in k(r), and hence a small change in the 
amplitude of the solution, but if s is of order a few then 
we may get a significant change in the phase <p(r). Indeed 
the phase difference can be written as 



A<p(r) 



Afc(r) dr + fc(r Biax )Ar n 
A[k 2 (r)} 



2k(r) 



■ dr. 



(A7) 



(The second term goes away because r max is a zero of k 2 .) 
The change in k 2 (r) has a contribution 2m e AE/h 2 if we 
change the energy, and another contribution — A[L(L + 
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l)]/r 2 if we change the angular momentum. Thus we 
have 



Aip(r) 



ml /2 AE f r 



dr 



y/2eh J r y/ 1 — 1 
hA[L(L + 1)] 



r" 1 

' max 



2y/2m e e 



dr 



r r^v r 



_=. (A8) 



It is easy to verify that for A[L(L + 1)] of order unity (it 
is 2 for S-P transitions and 4 for P-D transitions) and 
r/ao greater than a few, the second integral produces a 
phase shift of Acp(r) -C 1. Therefore we drop it. The 
function A<p(r) can be solved analytically but is most 
easily expressed through the following parametric form. 
Let us define the dimensionless function 



2r -3/2 



dr 



V^ 1 



(A9) 



so that A<p(r) = {ml' 2 rULAE / 2^/2 eh)r. Then by the 
substitution r = r max (l + cos?y)/2 we can derive r = 
?7 + sin 77, i.e. the relation between r and r is a cycloid 
function. Note that r = at r = r max and t = ir at 
r = 0. 

The matrix elements between two levels n 1 L and n' 1 P° 
depend on the integral 



K^Li r ) R n'- L P"(r)r dr; 



(A10) 



noting that L changes by 1 between the initial and final 
states, and that for small s = n! — n the normalizations 
of the wave functions are very similar, we may write 



Pn.n' . 



x r— cos <p(r) cosj^rj + A<p(r)] ^ ^ (An) 

JO V r max/ r - 1 

If we note that <^(r) is rapidly varying but A<p(r) is not, 
then the product of cosines can be averaged over several 
cycles to get 

R n .n>,L~- / r p=L dr. (A12) 

7Tr max J Q \Jr maSL /r - 1 

Changing variables to r gives 



2vr 



/ (1 + cos 77) cos dr. (A13) 
Jo 



The integrand is even in r so we may extend the range 
of integration down to —it and divide by 2. We may 
also replace the cosine by a complex exponential since 
the imaginary part is odd in r and hence vanishes. This 
gives 



(-iy 



-in 



(1 + cos rj)e inT dr, (Al4) 



where Q — 2 r^l 2 ^ AE / 2^/2 eTi. Now for large n, the 
energies are given by 



E 



2ao(n 



(A15) 



where Sl is the quantum defect for angular momentum 
L [52]. Note that Sl = for the hydrogenic case, but in 
helium there is a nonzero value due to the complicated 
physics occuring at small r (of order ao). The quantum 
defects for He i singlets are -0.1397 ( 1 5 < ), 0.0121 ^P ), 
and -0.0021 ( 1 L>) [H3, US]. Therefore we have AE w 
(e 2 /aon 3 )(s + Si — Sl) and r max ps 2aon 2 , which implies 
ft w s + Si - 5 L . Thus 



\S-1 



/ (l + cos^e^^-^^dr. 

J — 7T 



(-1) 

Rn,n',L ~ — ^ aon 

Ztt 

(A16) 

Thus we see that the radial matrix element is simply the 
Fourier transform of the cycloid function. This is consis- 
tent with semiclassical intuition since the cycloid is the 
classical trajectory of a particle in a Coulombic potential 
with very small angular momentum. The reduced ma- 
trix element required to compute Eq. lfl3|) is obtained by 
multiplying by the relevant angular factors: 

(n^Hdlln'V) « (-l) i >+ s - 1 Ll/ 2 e a n 2 / cyc ( S + ( 5 1L ), 

(A17) 

where we have introduced the shorthand S\l = Si — Sl 
and 

/cyc(O) = i- J (1 + cos77)e inr dr. (A18) 

The key result - the cancellation of contributions to 
AI27 for large n - comes from the following identity: 



E (-l)7cyc(* + £) 



(1 + cos 77) 



S— — OO 



dr 
27 



= £(1 + 003^-0(^ + 1)1 = 0, (A19) 

since 1 + cos 77 = when r is an odd multiple of ir. (Here 
IH is the sampling function.) 



APPENDIX B: STEADY-STATE LINE WITH 
FINITE LINEWIDTH 

In this Appendix we consider a simple analytic model 
for the steady-state line profile in the vicinity of a res- 
onance, i.e. an approximate solution to Eq. |4T|. This 
approximation is valid when the half-width of the "reso- 
nant" part of the spectrum, Ai/ cut , is small compared to 
the frequency diference between neighboring resonances 
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as well as compared to the thermal scale k^T r /h. As 
an example, for the He I 2 1 P°-1 1 S' line, we have used 
A^ cu t = 58THz; the frequency distance to the next al- 
lowed resonance (2> 1 P°-\ 1 S) is 452 THz; and the thermal 
scale is 110[(1 + z)/2000] THz. 

It is easily seen that the matrix element -M 27 for the 
i — > 1 : 5 two-photon process possesses a simple pole at 
each resonance v , v' = AE(n' L P°)/h. Therefore the two- 
photon decay rate, which is the square of the matrix ele- 
ment times phase space factors, can be written in a power 
series 



dT 



2", 



du 



E 

H=-2 



(Bl) 



where Av = v — v 1 i s _ n ii po . The power series cuts off at 
fi = — 2 because the square of a function with a simple 
pole can have a pole of no higher than the second order. 
It is easy to read off from Eqs. (|12p and (fl3]l that the 
leading term for the i — > n' 1 P° ^ l 1 S pole is 

a I/flO /lDo^ 3 /lDo 

l L S—n' L P° n' P°—%{-\ i k r/ \i 

*'- 2 = 108(2L + l)agtt° 

x (l 1 5||d||n /1 P }(n /1 P°||d||i)| 2 , (B2) 

where we have taken N(y\is-n n P°) *C 1 in the Wien 
tail of the CMB. Using the conversion from dipole matrix 
element to Einstein coefficient, and replacing the phase 
space density with its blackbody value, this can be re- 
written as 



1i,-2 



■^■i->n ll P°-^n n P'>->l 1 S 

4 7 r 2 (l-e~' uy "' 1 f°-* /fcBTl )' 



(B3) 



It follows from this and Eq. I|37p that the absorption cross 
section to level i is 



-f 7) = E 



where the leading order term is 

A-m' 1 P° ^ri' 1 P° ^ 1 1 S 



Qi,-2 = 



c 2 gi 



This could alternatively be written as 



(B4) 



(B5) 



32^^ 



-^n' 1 P°-l 1 S r n' 1 P°-*i- (B6) 



S-n 7l P° 



A similar argument shows that Eq. (|B6|1 applies to the 
resonance in the Raman scattering cross section corre- 
sponding to l 1 ^ ^ n' 1 ?" ^ i as well. 
Equation (|4l|) thus becomes 



dNo 
dv 



(B7) 



where A/L = Vp» /3%v-s an d 

K= ^ ' 



(B8) 



Expanding n as a power series, K = $^L_ 2 K ix^- V ^ '> we 
find that the lowest-order term is 



E^HCXii S 
77 W,-2 



(B9) 



In the final expression, the prefactor outside the sum is 
easily recognized as t$ /4ir 2 , where r$ is the Sobolev depth 
through the line. The sum is the total width of the n' P° 
level (which is the line width ru ne of n' 1 P°-l 1 S since 
the 1 X S level has negligible width) times the fraction of 
transitions from n' 1 P° that go to other excited states. 
Therefore the sum is rii ne /j nc and we may write 



K-2 



TsEline/i, 

4tt 2 



Av Yv 



(BIO) 



Note that the coefficient k goes to infinity on resonance. 
In principle this should be cut off by the Lorentzian width 
of the line (i.e. the pole displacement in ^27), and the 
resonance will also be widened by the Doppler width of 
the line. In practice as long as the line center is optically 
thick this subtelety does not matter: we will have A" = 
A/l at Av = 0. 

Our next objective is to solve Eq. i|B7|) for small Av. 
Here we take "small" to mean that we can work to first 
order in hAv/hsT? and the correction terms {«^}S=_i- 
We may begin by writing the solution, 



Ao = -X / A-^ALe-' 1 ^/^ du, 



(Bll) 



where X = cxp J ndv. The constant of integration in 
X is arbitrary (it trivially cancels out in obtaining Ao), 
while that of the integral in Eq. (|B11[) is determined by 
boundary conditions. We will separately solve for the 
Av > and Av < regions since X is singular at Av = 
0. The solution for X is 



/ Afi inc Av 
X = cxp h K-i m ■ 



Av cut 



E k " 

M=0 



A 4 



(B12) 

where the choice of denominator in the logarithm is 
arbitrary (but Av + is convenient). Substitution into 
Eq. <|B11|> gives 



e E~ ^a^+Vm 



x + £ Km a^J e - hA ^ k ^ dv. (B13) 
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Expanding this to first order in hAv/k^T? and {k ai }^ = _ 1 
gives 



X 1 + K-i In 



Av 

Av cut 



/j=0 



oo A u+l 00 

o""""^ A^ 



k u Av»+ 2 \ 

1 dv. (B14) 



The integral can be shown by direct differentiation to 
evaluate to 



a Ai/i ino /Ai/ 



, Av ^ A^ +1 
1 — K_iln— > k u 

cut h * , 



Az/ 



(B15) 



The photon phase space density on the red side of the 
line is easiest to obtain: since the term multiplying C in 
Eq. (|B16jl goes to infinity as Av — ► 0~, we must have 
C = 0. We thus have 



1 + ^A^ inC cA ^. nc/A „ cut Z' Ai/ Iinc \ 



V Afc'cut / 

(B17) 

Using the expansion of the exponential integral for small 
values of the argument, we find that if Avu ne <C A^ cut -C 
feTr/Zi, then 



where C is a constant of integration and E\ is the ex- 
ponential integral function. Therefore, to first order in 
K}£°=-1 and hAv/k B T r , 



fc B T r 1.78A^ 



(B18) 



^ ^A^We e -A!/ line // 



fe B T r 

I Pg-A>1iiie/Af 

x 1 + k_i In 



Si - 



AiAi, 



Av 



Av cnt 



(B16) 



where 1.78 = e 7 is the exponential function of Euler's 
constant. 
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